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We study dissipative inflation in the regime where the dissipative term takes a specific form, 
r = r((j!>), analyzing two models in the weak and strong dissipative regimes with a SUSY breaking 
potential. This system introduces three new parameters; two for the potential and one for the dis- 
sipative term. After developing intuition about the predictions from these models through analytic 
approximations, we compute the predicted cosmological observables through full numerical evolu- 
tion of the equations of motion, relating the mass scale and scale of dissipation to the characteristic 
amplitude and shape of the primordial power spectrum. We then use Markov Chain Monte Carlo 
techniques to constrain a subset of the models with cosmological data from the cosmic microwave 
background (WMAP three-year data) and large scale structure (SDSS Luminous Red Galaxy power 
spectrum) . We find that the primordial dissipative parameters are uncorrelated with the "late-time" 
cosmological parameters which describe the contents and expansion rate of the universe; the latter 
show no significant shift from the standard ACDM concordance cosmology and possess close to 
Gaussian posterior probability distributions. In contrast, the posterior distributions of the dissipa- 
tive parameters are highly non-Gaussian and their allowed ranges agree well with the expectations 
obtained using analytic approximations. In the weak regime, only the mass scale is tightly con- 
strained; conversely, in the strong regime, only the dissipative coefficient is tightly constrained. A 
lower limit is seen on the inffation scale: a sub-Planckian inflaton is disfavoured by the data. In both 
weak and strong regimes, we reconstruct the limits on the primordial power spectrum and show that 
these models prefer a red spectrum, with no significant running of the index. Despite having one 
extra degree of freedom in the fit compared to the standard ACDM model, the data does not display 
a preference for any of the dissipative models; the goodness-of-fit is comparable to the latter. We 
calculate the reheat temperature and show that the gravitino problem can be overcome with large 
dissipation, which in turn leads to large levels of non-Gaussianity: if dissipative inflation is to evade 
the gravitino problem, the predicted level of non-Gaussianity might be seen by the Planck satellite. 



I. INTRODUCTION 

In the original model of inflation, the early universe 
is hypothesized to undergo exponential expansion and 
supercool as a result of the radiation redshifting A 
subsequent, independent phase of reheating is required 
to return the universe to a radiation dominated state 

This original picture assumes that any interaction of 
the inflationary scalar field (called the inflaton) has a 
negligible consequence on infiation. It has, however, been 
shown that a scalar field interacting with a thermal bath 
leads to an additional friction term in the equation of 
motion Q. Hence, if the universe begins in a thermally 
excited state, the inflaton equation of motion should pick 
up a sizable friction term from thermal interactions 
The interaction results in radiation production, which 
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prevents the universe from supercooling. The significance 
of the extra dissipative term was independently realised 
in and later in p^ . where the inflaton dynamics were 
studied and the scenario was named warm inflation. 

Over the past several years, warm inflation has been 
shown to be a valid model of the early universe. The 
background dynamics of the dissipative inflaton has been 
modeled and the non-equilibrium thermodynam- 

ical problem has been studied extensively p^. [T^. [isl. [igI. 
[TtI . [isl . [l9| . The characteristic friction term has been 
calculated for sup ersymmetric models with a two stage 
decay process j20l[2]|. Cosmological implications of dissi- 
pative terms in h ybr id inflationary models have been ex- 
tensively studied [22l.[23j. Since the characteristic friction 
terms are not limited to finite temperature effects, the 
more general name dissipative inflation shall be adopted 
here. 

A study of the cosmological perturbations produced by 
dissipative infiation has identified interesting character- 
istics in the cosmic microwave background (CMB) power 
spectra from warm inflation (2^ ; a small level of running 
of the scalar spectral index is predicted, while tempera- 
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ture dependent coupling between the field perturbations 
can lead to oscillatory spectra. Another distinct char- 
acteristic of these models is t he p redicted level of non- 
Gaussianity, which is high [H, Ha [13, Hi • 

A naive best-fit parameter search has been previously 
completed using the first year WMAP CMB data for one 
model (2^ . demonstrating the existence of parameters 
which map this model of warm inflation to cosmological 
data. This paper aims to extend the previous work, by 
(a) studying several models of dissipative inflation, (b) 
using robust Markov Chain Monte Carlo (MCMC) tech- 
niques to constrain the model parameters in a subset of 
models, and (c) using two datasets which significantly 
improve upon the previous analysis, namely the WMAP 
three-year data [3C| and the SDSS Luminous Red Galaxy 
(LRG) power spectrum data [3ll |. 

This paper is organized as follows: Section [ll] intro- 
duces the dynamics of dissipative inflation, both for the 
background and perturbations. Section IIIII defines the 
amplitudes for both quantum and thermal fluctuations. 
Analytical approximations for the numerical models con- 
sidered are derived in Section IIVI split into three main 
regimes; weak (r <C 1, T < 7?), thermally weak (r <C 1, 
T > H) and strong dissipation (r » 1, T > iJ). The 
implications of analytic limits are discussed in Section 
fVl and the numerical simulation and cosmological con- 
straints are detailed in Section lVll The results are finally 
given in Section IVlIl 
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FIG. 1: The full SUSY breaking potential (solid line) used 
for numerical simulation and the m^(j}^ potential (dotted line) 
used in approximation. Here m = 0.95/^ = 0.95Mpi. 



p, can be written 

P 
P 



2-r ' yri ' ^Hr 
and the zero curvature Friedmann equation becomes 



II. DISSIPATIVE DYNAMICS 
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The dynamics of dissipative inflation are described in 
[2^, where a full discussion of the background and per- 
turbation solutions can be found. In an expanding, ho- 
mogeneous universe, the dissipative inflaton equation of 
motion is given by 



(3H + r)0 + 1/^ = 0, 



(1) 



where V{(f>, T) is the thermodynamic potential and 
T{(j},T) is the damping term due to interactions between 
the inflaton and surrounding flelds. Both the thermal 
correction to the inflaton potential and the damping force 
are directly due to the effect of the interactions. Issues 
regarding thermal corrections lead to some debate as to 
the viability of warm inflation [3^ . The corrections to the 
effective potential were calculated for one explicit model 
and were shown to be negligible [s^. Such corrections 
shall be ignored for the purpose of this study, in which 
we will consider a SUSY breaking potential of the form 



where Mp\ is the reduced Planck mass, 2.436 x lO^^GeV. 
Due to the damping term, the continuity equation ac- 
quires a source term, which indicates radiation produc- 
tion: 



Pr + 4:Hpr = r0 



(4) 



The temperature of the heat bath can be calculated using 
Pr = 7r^g*r^/30, where g*{T) is the effective particle 
number. The ratio of radiation production to expansion 
rate. 



= r/3H 



(5) 



is a useful parameter to distinguish between regimes of 
weak or strong dissipation. 

It is usual to assume a slowly-rolling inflaton, so that 
the second derivatives in the equations of motion may be 
ignored. In this limit, slow-roll is determined by a set of 
parameters: 
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The form of this potential is shown in Figure P]). 

Accounting for the radiation density, Pr, for a homoge- 
neous scalar field, the total density and pressure, p and 
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The first three are required to be less than 1 + r. 
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In the zero-shear, longitudinal gauge, 

ds^ ^ -(1 - 2ip)dt^ + a^{l + 2Lp)5^jdx'-dx^ , 
the perturbed Einstein equations lead to 



3if(y9+ (3if2 + fc^a"^) ^ 



(p - AHlp - (2H + LP 
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with 



6p — 4>S(j) + V,^5(l) + (j)'^ip + 6pr, 

5p = (pSij) — V,^S(t) + (j>'^ip + -6pr. 

o 

In addition, the perturbed equation of motion becomes 



T4,Lp - 2V^ip = 0. 



(8) 



Note that when T = r((/),T), the (p{Sr) term leads 
to a coupling of the perturbed radiation and inflaton 
fields. This non-trivial coupling is important in deter- 
mining the shape of the power spectrum of perturba- 
tions. Models with inverse-temperature dependent fric- 
tion terms have been shown to produce oscillations in the 
power spectrum, with diminishing amplitude for increas- 
ing wavenumber 2J| . While these models are phenomeno- 
logically interesting, we note that, in these models, each 
oscillation contains a minimum at Vk = 0, which result 
in zeroes in the power spectrum. Since these are not seen 
in the acoustic spectra, we assume they cannot reproduce 
observations. Further work is required to reduce the am- 
plitude of such oscillations in these models in order to 
convincingly match to data. In addition, friction terms 
with positive temperature dependence may not lead to 
oscillation. These models are beyond the scope of this 
paper, but is the topic of on-going research. In this work, 
we will restrict the form of the friction term to F = F((/)) . 
One can consider the gauge-invariant variable [H, [35l . 
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where v is the velocity perturbation. In the longitudi- 
nal (zero-shear) gauge, this variable is equivalent to the 
curvature perturbation in the comoving gauge [13] . Well 
outside the horizon, k <C ai7, 
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Typically, TZ quickly freezes outside of the horizon and 
the final amplitude can be approximated using horizon- 
crossing field values. The power spectrum for the curva- 
ture perturbation is defined by 

7'K(fc)(27rfc-i)3,5(3)(fc-fc') = {RkR-k') 



(11) 

(12) 



III. FLUCTUATION AMPLITUDES 

In supercooled inflation, fluctuations of the inflaton 
field are seeded by vacuum fluctuations inside the hori- 
zon. The amplitude is well known (38l . [s^ . [40| and given 

by 



sm = |. 



(13) 



In dissipative inflation when T > H, due to the ther- 
mal bath, thermal fluctuations in the fields are expected 
to dominate over their quantum counterparts. For the 
strong dissipative regime, r ^ 1, the amplitudes are cal- 
culated in [iB,[23|- III the next section, however, it will be 
important to consider cases where r « 1. The fluctuation 
amplitudes (with F consistently replaced by 3H{1 + r)) 
are given by 
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As in [24|, the sign of (50 in Eqn. (|T5|) has been chosen 
for consistency with cross-correlations. Notice that for 
the str ong regime, r ^ 1, these results return the val- 
ues in [2J|. In contrast with supercooled inflation, ther- 
mal fluctuations do not freeze at Hubble crossing, but at 
a time given by ka^^ = (FiJ)^/^ called thermal freeze- 
out [l^. When r > 1, the freezeout time always precedes 
the Hubble crossing time, at which ka~^ = H . We note 
that freezeout has been calculated more recently with an 
additional factor of \fij2 [23], but for the numerics we 
have not included this pre-factor and it does not affect 
our conclusions. 

It should be noted that, even for r < 1, thermal fluc- 
tuations may dominate over quantum ones, ii T > H. 

The fluctuations of the radiation field remains as in 
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Using the above terms for 5(f) in Eqn. (|12[) . we may 
calculate the spectral index, — 1 = dV-jz/dhik and 
running of this index, dus/dlnk [3, Slj, when F = F((/)). 
We require two second-order slow-roll parameters: 
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where the final equation is calculated at horizon crossing. 



For quantum fluctuations, using Eqn. (jl3[) . corrections 
due to the additional friction are found (r < 1): 

ris-l = -6e, + 27?* + n(8e*-27/, -2/3,) (18) 

4^ = 16e*?7* - 24e2 - 2^ + r, (646^ + 47, 
d In A; ' ^ 

-14e,/?* - 38?7,e* + 4/3,77, + 4C^) 
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For thermal fluctuations when r < 1, using Eqn. (jl4p : 

/21 3^ 3 
dris _,^2 21 3^ 3 3^ 



din A: 
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-ye,7?* + ^/3*r/* + J/3' + 3C' 
When r > 1, thermal fluctuations lead to: 
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Due to the constancy of the curvature perturbation out- 
side the horizon, the spectral index may be calculated at 
horizon crossing or freezeout, as indicated by an asterisk. 
The spectrum of gravitational waves is given by [4l| 
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If fluctuations in the radiation field exist, then isocurva- 
ture fluctuations should be expected. Their spectrum is 
given by [i^ : 



Pi 
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Note that the difference in prefactors with respect to [43 
arises from the definition of the scalar perturbation am- 
plitude. The scalar-to-tensor ratio and ratio of isocurva- 
ture perturbations can be defined by 
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which approximates closely to the full potential close 
to the minimum (see Fig [1]). A good fit occurs when 
m « 0.95/i. In addition, we shall restrict ourselves to the 
simple case in which 



r(0) = Fo 



(25) 



The slow-roll regime can be assumed, for which the 
equations of motion simplify: 
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and the slow-roll parameters are given by 
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For this potential, the index is given by: 

(4 — 6r) r < 1 (quantum) 



0.(27) 
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^r*) r < 1 (thermal) 
r > 1 (thermal) 



where we have assumed that (3 is negligible compared to 
e and 77 (which is really only true as (j) — > 0o)- Notice 
that for all three CcLSGS, Oil red index is expected. Due to 
the dominance of e (and 77) , we expect the running to be 
0((1 — ris)^), hence of the order of 10"^ or less. 



A. Background Values 

In the weakly dissipative limit, 1 — w 4e* leads to 

"8 Mpi 
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Conversely, in the strong limit. 
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For large Fq, C is small and, for all n, it is clear that 
either 



IV. ANALYTICAL APPROXIMATION 

In the next section, we will numerically evolve the full 
equations of motion, solving for both the background and 
perturbations. The full system is complicated and the 
results are non-intuitive. It is therefore useful to find an 
approximate analytic solution, which should provide a 
good guide for the expected parameter ranges. The full 
potential is the SUSY breaking one given in Eqn. ([2]). 
Near the minimum, it is possible to consider a shifted 
m^(/)^ potential: 



T/(0) = (00 
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We are interested in the case in which inflation occurs 
very close to 0o and hence we take this limit. Note that, 
in this regime, the friction is very large and hence a large 
number of e-foldings may be obtained even when the in- 
fiaton is close to its minimum. 

At the end of inflation, e « (1 + r). In fact, for strong 
dissipation, e sa 27- . Hence, in these two limits. 
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This agrees with the previous results. Eqn. ((28)) and ([30|) . 
where (j)f (pQ in both hmits. 

It is convenient to plot the dependence of r on Tq and 
TO, in order to determine the weak and strong regimes. 
For this potential, 
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This ratio at horizon-crossing and at the end of inflation, 
and r/, is dependent on 0* or respectively, which 
in turn depend on the parameters Fq, to, and as in 
Eqns (|29II32|) . Prom Eqns and we may take an 
average value of 0* w O.80o- For (pf we solve Eqn. ([32]) . 
Results for r* and r f for relevant ranges of Fq and to are 
plotted in Figure Q for n = 1. 

The radiation density can be calculated using the slow- 
roll approximation ^43||, 
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from which the temperature can be approximated. Sub- 
stituting rf into this equation leads us to the final reheat 
temperature, plotted in Figure ([3]) for n = 1. When 
n > 1 , the final values r f and Tf do not differ much from 
the n — I case, since (jjf w (pQ. 



From Eqn. 
ratio T/H: 



, it is also possible to approximate the 
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Note that when r <^ 1, T/H cx j-^/"* and the ratio is 
suppressed by r (although increases with r). When r 3> 
1, T/_ff cx r"^/"*, but since e w (1 — ns)r, the dependence 
on r cancels. 

Finally it is possible to estimate the number of e- 
foldings during infiation: 
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resulting in 
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In all the cases we shall now consider, the scalar-to- 
tensor ratio is set by the scale of the potential: 
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A very conservative upper limit is found when m/Mpi ~ 
10~^ and (po/Mpi ~ 10 (values considered later in the 
paper), such that 



Rg < 0.2. 



B. Weak Dissipation Amplitudes 

In the weak dissipative limit, the quantum fiuctuations 
dominate and the amplitude of fluctuations is given by 
Eqn. ((T3|). In this case, the amplitude of perturbations 
at horizon crossing is 
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Substituting « (1 — ns)(l + |?')/4, we find: 
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Note that the result is independent of n, since the am- 
plitude is independent of F (if we assume r is negligible). 
The range of values allowed given Vn and Us is shown 
in Figure g]). For typical WMAP values, to « 10"'^Afpi. 
We note that the scalar-to-tensor ratio, Rg, is calculated 
using Eqn. (PT|) and is given by the standard supercooled 
limit of Rg = 16e. The ratio of isocurvaturc perturba- 
tions is given by 

9 ^ 1/4 ^3/4^5/4 ^,^pj 



Ru 



2a 



6 X 10" 



507r 



\Mpi 



C. Thermally Weak Amplitudes 

In this regime, the friction term acts only weakly r < 1, 
but the thermal amplitude dominates over its quantum 
counterpart, such that the amplitude is given by: 
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For r < 1, from Eqn. (j33[) we require: 
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FIG. 2: The ratio of the friction terms at horizon crossing (r, , LHS plot) and end of inflation (r/, RHS plot), which determines 
weak (r < 1) and strong (r > 1) regimes. A value of <;/>« = 0.8<^o is assumed in Eqn. (|33} as detailed in the text. A value of 
n = 1 has been assumed and (po = Mpi (solid line). For reference, the line r, = 1 when (po ~ lOAfpi has also been plotted 
(dotted line). 
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FIG. 3: The final reheat temperature, Tf, over the full pa- 
rameter range considered. A value of n = 1 has been assumed 
and (fio = Mpi. 



and imposing T/H > 1 (with r < 1), wc find 

n 3 
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These constraints on the parameters m and Fq are shown 
in Fig. m when (j>o/Mp\ = 1, 10. Also superimposed are 
the required values of (/)o which satisfy Eqn. (|43|1 for typ- 
ical WMAP values. The ratio of isocurvaturc perturba- 
tions is given in this regime by 

3r e 



D. Strong Dissipation Amplitudes 

For the perturbations in this regime, the new friction 
term dominates (r ^ 1) and, assuming T > H, the am- 
plitude is given by: 
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where we have used pr = 7r^p,T^/30 = aT^ and (1 — 
TT-s) « 3e*/4r'*. There is still some field dependence in 
this equation, due to the form of r(0), but from Eqn. 
pop wc sec that this dependence is very weak and can be 
ignored. For this reason, the dependence on the index n 
is also very weak. 

It is immediately obvious that, in this strong limit, the 
amplitude in Eqn. (|46p is independent of the mass m. In 
fact, 
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where the last equation has used — 100 (a sa 32.9). 
Figure ([5]) shows how Fq depends on realistic values of 
Vn and when n — \. For typical WMAP amplitudes, 
Fo ~ lO^^Mpi. The ratio of isocurvaturc perturbations 
in the strong regime is given by (4^ : 

1 e 
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FIG. 4: Expected parameter values in the weak regime (r <^ 1) when n — 1. When T < H (LHS plot), the scale of the 
potential, m, is constant, given in terms of the amplitude V-jz, and index ris as in Eqn. (|42|) . When T > H (RHS plot), the 
three parameters are related by Eqn. (|43|l and are plotted for typical WMAP amplitudes. In addition, the constraints for r < 1 
and T > H are also shown for <j)o — A/pi (solid) and (fio ~ lOAfpi (dotted). Note that when T > H, a, further constraint comes 
from iV > 60 (or 0o > Afri). 
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FIG. 5: Expected parameter values in the strong regime 
(r ^ 1) when n = 1. In this strong regime. To is constant 
and is given in terms of the amplitude Vn, and index Ua by 
Eqn. (HZll. 



V. IMPLICATIONS FOR PARTICLE PHYSICS 

The analytic results lend some insight into the allowed 
physical scales of the inflationary potential. In the weak 
dissipation regime, as expected for supercooled inflation, 
the amplitude of the perturbations, as seen by observa- 
tion, merely sets the scale of the potential for the inflaton 
field, as seen in Eqn. (|42)) . Conversely, however, no limits 
are placed on the scale of the inflaton in the strongly dis- 
sipative regime, whilst the amplitude of the friction term 
is fixed. 



One of the most stringent constraints on these models 
will be the reheat temperature, which is of order 10^*^- 
10^"* GeV in these models. Counter-intuitively, weak dis- 
sipation normally leads to a larger reheat temperature. 
In order to evade an overproduction of gravitinos, the 
reheat temperature should be less than 10^ GeV Hj, 
(for a review see [1^). From Figures ([2]) and it is 
clear that this requires r ^ 10'^, which in turn requires 
Fo/A/pi > 10"^ for weak dissipation or m/Mpi < 10"^° 
for strong dissipation. Alternatively, one must consider 
a mechanism to dilute the gravitinos after infiation ends, 
such as late-time entropy production. 

Additionally, due to high levels of predicted non- 
Gaussianity, extra constraints are placed on the param- 
eters. From prj . the non-linearity parameter, /^^'' (m 
the equilateral limit) is given by 

/^r--151n(l + ^)-| (48) 

Present constraints on the non-linearity parameter give 
-256 < f^l'^ < 332 at 95% ^ and result in r < 10^". 
This places a rather loose lower bound on m, 0o and 
Fq, as given in Eqn. (|33p . However, this constraint is 
expected to improve with data from the Planck satellite; 
a 1-a limit on non-Gaussianity |/]^^''| < 66.9 is forecast 
for the latter [Hi. This corresponds to r « 1000 [H]. If 
the dissipative regime evades the gravitino problem, the 
level of non-Gaussianity should be seen by Planck. 

VI. NUMERICAL SIMULATION 

Now we will discuss the numerical simulation which 
forms the basis of our cosmological constraints. It is im- 
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portant to note that we will only use the previously de- 
rived analytic results to cross-check the numerical code. 
No analytic or slow roll approximations are used in the 
numerical computations of the cosmological observables 
which are compared with the data. 

It is convenient to adopt a general form for the friction 
term 

ri^) = ro(^-^y . (49) 

In order to check the above analytical approximations 
and identify the effect of higher powers of (f>, we will con- 
sider two numerical cases: 

• Case I: n = 1 

• Case II: n = 2. 

The case in which ?i = (a constant friction term) is 
unphysical and will be ignored. 

We choose to study both strong and weak dissipative 
regimes for both cases I and II. Of course, in the weak 
regime, there are two separate scenarios: T > H and 
T < H ("thermally weak" or "weak"). We therefore 
study a total of six models: weak, thermally weak and 
strong, for both n= 1,2. 

We expect the models to give a flat spectrum of cur- 
vature perturbations, with a small, inherent amount of 
blue to red running [2^. Although the running is 
small, we are interested in how this running leads to a 
constraint of the parameters. The chosen inflaton poten- 
tial originates from soft supersymmetry breaking and is 
given by Eqn. ([2]). Note that slow roll supercooled infla- 
tion cannot occur with this potential, since the standard 
slow-roll conditions are difficult to satisfy. 

A. Dissipative Inflation Code 

We consider the evolution of perturbations for both 
cases I and II in weak, thermally weak and strong 
dissipative scenarios. The background dynamics fol- 
low Eqns. mil), while the perturbations are given by 
Eqns. dZKHl). The code follows Ref. ^ with a few mod- 
ifications, detailed here. 

For T < H, the perturbations freeze at the Hubble 
radius {k = aH) with an amplitude given by Eqn. (jl3p . 
When T > H, thermal fluctuations dominate; freezeout 
occurs earlier at a time k = a{THy^^ and the ampli- 
tudes are given by Eqns. (|14j|16p . Wavenumbers k/a are 
normalised by the Hubble length at the end of inflation: 

■"(''^)-'°(-''-")— <») 

where T14 x 10^^ GeV is the temperature at the end of 
inflation. 

We separate the weak and strong dissipative scenario 
from the weakly thermal regime. This latter regime leads 



to a tightly constrained region of allowed parameter space 
and will be discussed separately in the next section. 

For what follows, we consider the two separate regimes 
of weak and strong dissipation, r < 1 and r > 1, with 
quantum and thermal fluctuations respectively. From the 
analytical constraints in Section pvp . we may assume 
numerical priors on the parameters, /i and Fq: 

r - 6.1 < logio (M/il/pi) < -5.7 
[ -10.1 < logio (Fo/A/pi) < -5.7 

r -10.4 < logio (mA^pi) < -7.1 
r » 1 <^ 

[ - 7.5 < logio (Fo/il/pi) < -5.7 

Analytically, there is no constraint on 0o from ei- 
ther strong or weak dissipation. However, from Sec- 
tion (|IV Ap . it is clear that 0o sets the total number of 
e-foldings. For numerical convenience, we impose a prior 
on 00 using the number of e-foldings: we impose an upper 
limit of < 10^, which requires log^o (</'o/-^Pi) ^ 2 for 
both cases I and II (also both strong and weak), which 
is greatly above normally considered values. 

There is a heuristic reason for why A^ < 10^ leads to 
the same cut-off in both strong and weak cases. In all 
the cases considered, the friction term increases mono- 
tonically with cj). Therefore, at the top of the potential 
(the beginning of inflation) , F is small and we always have 
a weak dissipative regime. From Eqn. p7p . it may be 
observed that, at this time, iV « (0o/Mpi)^. Of course, 
as r increases, the friction term becomes more important 
earlier and earlier and N will depend on Fq and m as 
well as 00. 

For strong dissipation, imposing a cut-off on A^ also 
imposes a prior on m as seen in Eqn. (|37p . This lower 
bound will be seen in the results later. 

A natural lower prior on 0o is imposed due to the 
requirement that cosmological modes cross the horizon 
during inflation (i.e. N > 60). 

As established in Section [TVl the isocurvature and ten- 
sor amplitudes for the ranges we consider are highly sub- 
dominant to the scalar amplitude and are negligible con- 
sidering the quality of data we compare to. Therefore we 
do not compute them numerically. 

B. Curvature Perturbation across the Horizon 

The curvature perturbation, TZ, which defines the am- 
plitude of the power spectrum, was shown to be approx- 
imated by H5(j)/(p^. Although it is usually assumed that 
the curvature perturbation does not deviate from its hori- 
zon crossing value, there is some time-dependence. Due 
to the additional terms in Eqn. ([8]), Scf) can oscillate, shift- 
ing the final value of the curvature perturbation. Hence 
the final value, TZf is not always closely approximated by 
7^*. 
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71 = 1 


n = 2 


A 


5.7 


176.7 


B 


3.24 


3.23 


C 


-0.19 


-0.69 


V 


-1.74 


-1.73 



TABLE I: Coefficients for the thermally weak cases. The 
empirical form of the amplitude is given in Eqn. (|51|l . Models 
are shown with n= 1,2. Note that the amplitude decreases 
when n increases due to the inverse dependence on F. The 
shift in C is due to the weak dependence on (j>Q. 



When r ^ 1, Eqn. ^ reduces to the supercooled case 
and the curvature perturbation does not change signifi- 
cantly across the horizon as expected. When r ^ 1, the 
friction acts to damp the oscillations and again no sig- 
nificant change in is seen. However, a slight shift in 
the above analytical estimates for m and Tq is expected 
numerically. 

In the thermally weak regime, however, the final curva- 
ture perturbation (as calculated numerically) is not well 
approximated by the freeze-out values. In this thermally 
weak regime, the final amplitude is not given by the 
expression in Eqn. although wc have numerically 

checked that the values at freeze-out are given by this 
expression. It is not possible to approximate the time- 
dependence of S(j) and instead we empirically fit the final 
amplitudes to an assumed form. Numerically, we predict 
a form 



I m \ I 



\Mp\J \M] 



PI 



/ i 



PI 



(51) 



where the coefficients are given in Table U for n = 1,2. 
Due to the small numerical range of (po considered nu- 
merically, < \ogii^{(t)[) / Mp\) < 2, the exact depen- 
dence on (j)Q is difficult to quantify. However, note that 
Vn^^^ cx {vn? /Tq). Therefore, for a given value of 00, the 
region of allowed parameter space in the log(TO)-log(ro) 
plane lies on the straight line, with gradient close to 2. 

Numerically, we also find that, close to this line, the 
amplitude can be predicted using Eq. (|5ip . but that the 
spectral index varies greatly. We therefore find that the 
width of this line of fit (as shown in Fig. (5]) is small, 
although it increases slightly with n, leading to a very 
tightly constrained fit. The best likelihood values we 
found are around 21n£ w 5392, which (as we shall see) 
are significantly worse than ACDM and the remaining 
cases we study. In addition, the reheat temperature in 
these models lies around 10^'^ GeV. Due to the narrow- 
ness of the constraint, wc will not consider this case in 
the following cosmological parameter estimation. 



C. Cosmological Parameter Estimation 

The fundamental observables of the dissipativc infla- 
tion scenario are the same as in standard ACDM cosmol- 



o 




-9.5 -9.0 -8.5 -8.0 
log,o(m/Mp,) 



-7.5 



FIG. 6: The line of fit for the thermally weak regime, as found 
numerically for n = 1 (solid line) and n = 2 (dotted line) . An 
empirical perturbation amplitude, Vn, is found as given in 
Eqn. (|5ip with coefficients in Table H] Note that the gradient 
of this fit is 2 as discussed in the text. These lines do not 
change for varying values of 0o within numerically considered 
ranges. 



ogy. These include the CMB angular power spectrum 
and the galaxy power spectrum. We modified the Boltz- 
mann code CAMB [4^ to calculate these observables; 
once the initial conditions are specified, the power spec- 
trum observables are computed from a separate primor- 
dial power spectrum routine as described above, which 
bypasses the usual CAMB code without significantly in- 
creasing the computational time. Using this code, we 
performed cosmological parameter estimation for the four 
scenarios specified by weak and strong cases I and II. 

We use a Markov Chain Monte Carlo (MCMC) tech- 
nique [13, [m, [H, m, m, [Sll to evaluate the likeli- 
hood function of model parameters. The MCMC is used 
to simulate observations from the posterior distribution 
V{0\x), of a set of parameters 6 given event a;, obtained 
via Bayes' Theorem, 



V{d\x) 



v{x\e)v{e) 
jr{x\e)vie)de' 



(52) 



where P{x\6) is the likelihood of event x given the model 
parameters 6 and V{6) is the prior probability density. 
The MCMC generates random draws (i.e. simulations) 
from the posterior distribution that are a "fair" sample 
of the likelihood surface. From this sample, we can es- 
timate all of the quantities of interest about the poste- 
rior distribution (mean, variance, confidence levels). A 
properly derived and implemented MCMC draws from 
the joint posterior density V{0\x) once it has converged 
to the stationary distribution. We use eight chains per 
model and a conservative Gelman- Rubin convergence cri- 
terion (s^l to determine when the chains have converged 
to the stationary distribution. 
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For our application, denotes a set of cosmolog- 
ical parameters. We then use a modified version 
of the CosmoMC code [s^ to determine constraints 
placed on this parameter space by the WMAP three- 
year cosmic microwave background data [sot and 
the SDSS Luminous Red Galaxy (LRG) galaxy 
power spectrum data (3l| . In our analysis, we 
take the parameter set {loi, = Q,bh?, Wm = ^■m.h'^, 
Oa, t, logio (-^o/Mpi) , logio W^Pi) , logio (ro/Mpi)}. 
Here, 9a is the angular size of the acoustic horizon and 
functions as a proxy for the Hubble constant Hq = lOO/i 
km/s/Mpc or ri„i. The universe is assumed to be 
spatially flat. Constant priors are assumed over the 
previously specified parameter set, subject to the mi- 
crophysical cuts applied in the dissipative inflation code 
as described above. Wc marginalize analytically over 
the linear bias factor b and the non-linearity parameter 
Qn\ of the SDSS LRG data as is done normally in the 
CosmoMC code. 



VII. RESULTS 

In this Section, we highlight and discuss what we con- 
sider to be the most important results of our cosmological 
analysis. Further details can be found in Figures ([HHSl) 
and their captions, and in Tables [H] and IIIII It should 
be noted that wc did not consider the thermally weak 
case for cosmological parameter estimation, due to the 
narrowness of the constraint. 

Table HIl compares the goodness of fit of the dissipative 
cases (weakly, thermally- weakly and strongly dissipative) 
studied with the "minimal" 6-parameter ACDM model. 
Since there is no improvement in the best fit even 
though the dissipative models have an extra degree of 
freedom, there is no indication that the data prefers the 
dissipative model over the ACDM case. 

Table Hill shows the ID constraints on the "late-time" 
cosmological parameters for the weak and strong cases 
after all other parameters have been marginalized over. 
The marginalized probability distribution functions for 
these parameters are very close to Gaussian, and they do 
not shift from the standard ACDM concordance values; 
they are all minimally correlated with the "early-time" 
dissipative parameters. We also expect this to generalise 
to the thermally weak scenario. 

Figures (O and ([8]) show the joint 2D constraints on 
the "primordial" dissipative parameters for the weak and 
strong cases respectively. Unlike the "late-time" cosmo- 
logical parameters, these parameters have highly non- 
Gaussian likelihood surfaces and are also very correlated. 

As predicted by the analytics in Section IIVBI for the 
weak cases, only fi is constrained; both Fq and (f>o ex- 
tend past their numeric priors and there is no correla- 
tion between these parameters. For both cases I and II, 
log^o {fJ-/Mpi) « —5.27, close to the analytic prediction. 
Additionally, the results are almost exactly independent 
of n (i.e. case independent) as expected. The sharp bend 



in the contour (for large Fq) occurs due to the increase in 
the friction ratio, r, at these values (seen on the LHS of 
Figure ([2])), as shown by Eqn. and is hence slightly 

dependent on the form of F in this region. 

There is also a noticeable shift in the mass, m, when 
(f)Q ^ Mp\. The low (f>Q values occur when the total num- 
ber of e-foldings, N, is just large enough (i.e. ~ 60), 
as seen by Eqn. (|37p . For these values, horizon crossing 
occurs very close to the top of the potential. Addition- 
ally, when (j)Q — > Mpi from above, the slow-roll factor, e, 
increases and therefore the spectral index decreases. In 
order to match observation, this shift must be therefore 
balanced by a shift in m. The case in which cosmologi- 
cal scales cross the horizon at the top of the potential is 
unphysical and is not of relevance. 

For strong dissipation, as predicted in Section IIVDI 
only Fq is tightly constrained. When r ^ 1, 0* w 0o, 
and logio (roA/pi) ~ —5.55 as given by Eqn. (|T7)) with 
almost no (j)o dependence. This may be compared with 
the numerical results of [1^ , in which a reasonable fit was 
found for [m/Mpi,(j)o/Mpi,ro/Mpi] « [10-^ 0.8, 10"^] 
(compared to WMAP I data). As r ^ 1 (^ increasing), 
the correlation between Fq and 00 is non-negligible and 
can be seen numerically. There is an upper constraint on 
fi due to the assumption of r > 1, but the lower limit 
should extend ad infinitum; the lower bound is derived 
from the N < 10^ bound we imposed. The upper bound 
of (j)o/Mp\ < 2 is also due to this e-folding constraint, as 
predicted. 

It should be noted that for both weak and strong cases, 
there is a lower bound of log^o (</'o/-^Pi) > 0, showing 
that the data disfavours 00 < Mp\. The constraint is 
slightly stronger for the weak case. This is due to the 
increase in the slow-roll parameter, e, which leads to an 
increase in the tilt (the index becomes redder) . As this in- 
dex become redder than the data, the fit becomes worse. 

Figure © shows the primordial scalar power spectra 
P{k) reconstructed from the MCMC runs using WMAP 
3 year and SDSS LRG data for dissipative inflation. For 
comparison, the reconstructed scalar power spectrum re- 
sults of Figure (10) of Ref. [l^l (see bottom right panel) 
are shown. The latter reconstruction was done using 
WMAP 3 year data and the SDSS main galaxy sample 
power spectrum, under the assumption that the primor- 
dial fluctuations are seeded by the standard single-field 
slow roll inflation mechanism that additionally satisfied 
a minimal "sufficient e-folds" requirement that solves the 
cosmological flatness and horizon problems. The dissipa- 
tive constraints are tighter than the single-field slow roll 
constraints for two reasons: (a) the SDSS LRG sample 
is more constraining than the SDSS main galaxy sam- 
ple, and (b) the single-field slow roll analysis marginal- 
izes over the shape of the inflationary potential using a 
Hamilton- Jacobi formalism [5^. 

Similar to the single-field slow roll constraints of ISTl . 
and consistent with the empirical power-law fits of [30| 
and [3l|, the dissipative constraints on P{k) show a dis- 
tinct preference for a red power spectrum. However, 
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Model 


Best fit -21n£,nax (WMAP+SDSS LRG) 


ACDM 


5374.58 


Weak I 


5374.78 


Weak 11 


5374.75 


Thermally Weak 


~5392 


Strong I 


5375.19 


Strong II 


5376.45 



TABLE II: The best fit chi square, defined as ~ —2 In £max 
(where >Cmax is the maximum likelihood with respect to the 
WMAP 3 year data and the SDSS LRG galaxy power spec- 
trum data) for the standard ACDM model and the four dissi- 
pative inflation scenarios analyzed in this work. The ACDM 
model gives a slightly better (lower) for this dataset than 
any of the dissipative inflation scenarios considered, while con- 
taining one less degree of freedom, and hence the data does 
not exhibit a preference for any of the dissipative models. 



there is a systematic difference in the reconstructed shape 
compared to the single-field slow roll case. The red spec- 
tra are predicted by the analytics in Eqn. ([28]) for both 
weak and strong cases. Due to the r-suppression in the 
strong cases, the tilt is expected to be less red than for 
the weak cases and this is what we see in the reconstruc- 
tion. No significant running is seen as expected. 

Numerically we have checked that the reheat temper- 
atures from the simulated models and the analytical ap- 
proximation agreed well. For example, the strong chains 
have reheat temperatures around Tfj « 10^^-lO^^GeV. 

Since the numerical code docs not simulate supercooled 
inflation, the constraints in this paper do not display a 
supercooled limit. It should be noted, however, that it 
is impossible with this potential to obtain enough su- 
percooled inflation to solve the standard problems (hori- 
zon, flatness); the parameters never reach a "supercooled 
limit" in the MCMC. 



VIII. CONCLUSION 

We have studied the dissipative inflationary regime 
when the dissipative term takes a specific form, F = F(0). 
Two forms of F have been analyzed in the weak, ther- 
mally weak and strong dissipative regimes, assuming a 
SUSY breaking form for the potential. Our system in- 
troduces three new parameters; two for the potential and 
one for the dissipative term. 

We derive analytic formulae which relate the mass scale 
and scale of dissipation to the characteristic amplitude of 
the primordial power spectrum and its shape. In the 
weak regime, the level of dissipation is unconstrained 
but the mass scale is set around 10~^Afpi (using typi- 
cal WMAP values). Conversely, in the strong regime, 
the mass scale is unconstrained, but the dissipative co- 
efficient is tightly constrained around 10~^Mpi. In the 
thermally weak regime, a non-trivial relation between all 



three parameters is found for the freeze-out value of the 
perturbation amplitude. Oscillation in Scj) across the hori- 
zon, in this latter case, means that the final amplitude 
cannot be analytically approximated and an empirical fit 
has been found. 

Using a numerical simulation of dissipative inflation 
to calculate the observable primordial power spectrum 
exactly, we explore constraints on the dissipative param- 
eter space given by the WMAP 3 year data and the SDSS 
LRG power spectrum using Markov Chain Monte Carlo 
techniques. We find that the values of the cosmologi- 
cal parameters do not shift significantly from the stan- 
dard ACDM concordance values. The constraints on the 
dissipative parameters show large correlations and agree 
well with expectations from the analytic understanding; 
the mass scale alone is constrained for weak dissipa- 
tion and only the strength of dissipation is constrained 
in the strong dissipative regime. In both cases, a sub- 
Planckian infiaton is disfavoured by the data. In both 
weak and strong regimes, we reconstruct the primordial 
power spectrum and show that these models prefer a red 
spectrum, with no obvious running of the spectral index. 
Despite the addition of an extra degree of freedom com- 
pared to the minimal concordance cosmology, the data 
does not display a preference for any of the dissipative 
models; the best fit is comparable to the standard ACDM 
model. 

In the models considered in this paper, the reheat tem- 
perature is between Tr « 10^°-10^^ GeV, which does not 
satisfy present constraints for gravitino production. This 
temperature can be reduced by raising the ratio of dis- 
sipation to the Hubble parameter. We also comment on 
the level of non-Gaussianity predicted in these models 
and relate them to our dissipative parametrization. 

This analysis does not consider a thermal dissipative 
term (F 7^ F(T)), which is expected to lead to oscilla- 
tions in the primordial power spectrum. These terms are 
beyond the scope of this work and will be the focus of 
future work. 
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Parameter 


Weak I 


Weak II 


Strong I 


Strong II 




U.UZZO'l^Q QQQ44 


022S4+°-<"""'5 


0.02302 ± 0.00047 


02304+" ""''** 




0.1066 ±0.0042 


0.1066 ±0.0043 


0.1070 ±0.0044 


0.1070 ±0.0044 


r 


0.092 ± 0.026 


0.091 ± 0.026 


0.109 ±0.025 


11 0+°°^^ 

U.iiU_o.028 


h 


U. ( Z»_Q QJ^g 


U. ( Z»_Q QJ^g 


726+" "^^ 

'^"-0.016 


742+f'0i6 

'^^-0.017 



TABLE III: Constraints on the "late-time" cosmological parameters from the WMAP and SDSS LRG data-sets (mean, upper 
and lower 68% CL, marginalizing over all other parameters), for the four dissipative inflationary scenarios described in the text. 
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FIG. 7: Joint 2D 68% (dark) and 95% (light) CL constraints on the "primordial" parameters of the weak dissipation case I 
(left) and case II (right). 
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[srj . The range of scales spanned by WMAP and SDSS LRG data (which were used in the fit) and the smaller scales covered 
by Lyman-Q data (which were not) are shown for reference. See text for discussion. 



[36] J. M. Bardeen, Phys. Rev. D22, 1882 (1980). 

[37] D. H. Lyth, Phys. Rev. D31, 1792 (1985). 

[38] A. Vilenkin and L. H. Ford, Phys. Rev. D26, 1231 (1982). 

[39] A. D. Linde, Phys. Lett. B116, 335 (1982). 

[40] A. A. Starobinsky, Phys. Lett. B117, 175 (1982). 

[41] A. R. Liddle and D. H. Lyth, Phys. Lett. B291, 391 

(1992), astro- ph/9208007. 
[42] A. N. Taylor and A. Berera, Phys. Rev. D62, 083517 

(2000) , astro- ph/0006077. 

[43] A. N. Taylor and A. R. Liddle, Phys. Rev. D64, 023513 

(2001) , astro-ph/0011365. 

[44] K. A. Olive, Phys. Rept. 190, 307 (1990). 
[45] M. Y. Khlopov and A. D. Linde, Phys. Lett. B138, 265 
(1984). 

[46] P. Binetruy, Super symmetry: Theory, experiment and 

cosmology (Oxford Univ. Pr., 2006). 
[47] P. Creminelli, L. Senatore, M. Zaldarriaga, and 

M. Tegmark, JCAP 0703, 005 (2007), astro-ph/0610600. 
[48] K. M. Smith and M. Zaldarriaga (2006), astro- 

ph/0612571. 



[49] A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 

538, 473 (2000), astro-ph/9911177. 
[50] N. Christensen and R. Meyer (2000), astro-ph/0006401. 
[51] N. Christensen, R. Meyer, L. Knox, and B. Luey, Class. 

Quant. Grav. 18, 2677 (2001), astro-ph/0103134. 
[52] L. Knox, N. Christensen, and C. Skordis, Astrophys. J. 

563, L95 (2001), astro-ph/0109232. 
[53] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002), 

astro-ph/0205436. 
[54] A. Kosowsky, M. Milosavljevic, and R. Jimenez, Phys. 

Rev. D66, 063007 (2002), astro-ph/0206014. 
[55] L. Verde et al., Astrophys. J. Suppl. 148, 195 (2003), 

astro-ph/0302218. 
[56] A. Gelman and D. Rubin, Statistical Science 7, 452 

(1992). 

[57] H. Peiris and R. Easther, JCAP 0610, 017 (2006), astro- 
ph/0609003. 

[58] H. Peiris and R. Easther, JCAP 0607, 002 (2006), astro- 
ph/0603587. 



